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ABSTRACT 

We use a comprehensive database of black hole masses (Mbh) and nuclear luminosities to investigate the re- 
lationship between radio emission and Mbh- Our sample covers a wide range of nuclear activity, from nearby 
inactive nuclei to classical Seyfert 1 nuclei and luminous quasars. Contrary to some previous studies, we find 
that the radio continuum power, either integrated for the entire galaxy or isolated for the core, correlates poorly 
with Mbh- The degree of nuclear radio loudness, parameterized by the radio-to-optical luminosity ratio R, also 
shows no clear dependence on Mbh- Radio-loud nuclei exist in galaxies with a wide range of Mbh, from ^10^ 
Mq to a few X 10^ Mq, and in a variety of hosts, from disk-dominated spirals to giant ellipticals. We demonstrate 
that R is strongly inversely correlated with L/L-e, the ratio of nuclear luminosity to the Eddington luminosity, and 
hence with mass accretion rate. Most or all of the weakly active nuclei in nearby galaxies are radio-loud, highly 
sub-Eddington systems that are plausibly experiencing advection-dominated accretion. 

Subject headings: black hole physics — galaxies: active — galaxies: nuclei — galaxies: quasars — galaxies: 
Seyfert — radio continuum: galaxies 



1. INTRODUCTION 

Efforts to search for massive black holes (BHs) in the cen- 
ters of nearby galaxies have made rapid progress in the last few 
,; years (Kormendy & Richstone 1995; Richstone et al. 1998; 
• van der Marel 1999; Ho 1999a; Kormendy & Gebhardt 2001). 
Kinematic observations of a significant number of galaxies 
have yielded evidence for central dark objects with masses 
- 10''-10''Mo, which can be plausibly interpreted as massive 
BHs. In two cases, namely the center of the Milky Way and 
NGC 4258, the extraordinarily high density of the dark mat- 
ter appears to rule out all reasonable astrophysical alternatives 
to a single collapsed object (Maoz 1998). To date, BH masses 
are available for ^40 inactive or weakly active galaxies from 
observations of spatially resolved kinematics and for a com- 
parable number of bright active galactic nuclei (AGNs) from 
"reverberation mapping" of their broad-line regions (see § 2.1 
and Table 1), sufficient to encourage preliminary examinations 
of statistical relationships between BH masses and global prop- 
erties of the galaxies. Two correlations have emerged for the 
weakly active galaxies: (1) Mbh correlates with Lbuise, the opti- 
cal luminosity (cx mass) of the bulge component of the galaxy 
(Kormendy 1993; Kormendy & Richstone 1995; Magorrian et 
al. 1998; Ho 1999a; Kormendy & Gebhardt 2001), and (2) Mbh 
correlates even more strongly with a^, the luminosity-weighted 
stellar velocity dispersion within the bulge effective radius 
(Gebhardt et al. 2000b; Ferrarese & Merritt 2000). By contrast, 
Mbh depends little on disk properties (Kormendy et al. 2001). 
These empirical relations hold great promise for furthering our 
understanding of the formation of massive BHs, the formation 
of galaxy spheroids, and the apparent close connection between 
the two. 

Another correlation which has received considerable atten- 
tion is that between Mbh and radio emission. Franceschini, 
Vercellone, & Fabian (1998) compiled a sample of the dozen 
BH masses known as of mid- 1997 and showed that they evi- 
dently scale with radio luminosity, approximately of the form 
Li-ad oc Mg^, where Lrad is measured at 5 GHz (6 cm). Sur- 



prisingly, the correlation is as tight, if not tighter, for the to- 
tal radio emission integrated over galactic scales than it is for 
the core emission alone. They (see also Di Matteo, Carilli, & 
Fabian 2001) argue that the functional form of the Liad-MflH 
relation arises naturally if the BH accretion takes the form of an 
advection-dominated accretion flow (ADAF; see Narayan, Ma- 
hadevan, & Quataert 1998b and Quataert 2001 for reviews), fu- 
eled by hot gas from a large-scale spherical inflow. Moreover, 
the steepness of the relation gives it considerable leverage in 
predicting Mbh efficiently from L^d, a readily available quan- 
tity. Several subsequent studies similarly considered related 
samples of nearby galaxies, but their results have been mixed 
(Yi & Bough 1999; Salucci et al. 1999; Laor 2000; Di Matteo 
et al. 2001). Discussion of the apparent dependence of Mbh on 
radio luminosity and the radio-loudness parameter {R) recently 
has also surfaced in the context of more active galaxies such as 
Seyferts and quasars (McLure et al. 1999; Nelson 2000; Laor 
2000; McLure & Dunlop 2001; Lacy et al. 2001; Gu, Cao, & 
Jiang 2001). 

This paper reexamines the relationship between radio emis- 
sion and Mbh in hght of the most up-to-date samples of ac- 
tive and weakly active galaxies with reliable BH mass measure- 
ments. We show that the loose trend between integrated radio 
luminosity and Mbh is largely indirect, a consequence of more 
fundamental correlations between radio luminosity and bulge 
mass on the one hand, and between bulge mass and Mbh on the 
other. The distribution of nuclear radio luminosity as a func- 
tion of Mbh is more physically grounded, but its large scatter 
renders it ineffective as a statistical tool to predict Mbh- The R- 
Mbh relation disappears altogether when one considers AGNs 
with a broad range of intrinsic luminosities. However, we find 
that R is strongly related to the mass accretion rate. 

2. THE DATABASE 

We begin with a fairly detailed documentation of the data 
used in the subsequent analysis. Particular attention is paid to 
the source of the Mbh measurements (Table 1) and photometric 
parameters (Table 2). 



2 



HO 



TABLE 1: Galaxies with Black Hole Masses 



Galaxy 


Hubble 


T 


Spectral 


cz 


D 


Ref. 






Method 


Ref. 


Name 


Type 




Class 


(km s^^) 


(Mpc) 




(M 


a) 








(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 




(9) 


(10) 


3C 120 (Mrk 1506) 


SO: 


-1.0 


SI 


9896 


137.8 


1 


2.3 


X 


10^ 


R 


1 


3C 390.3 (VII w 838) 


E? 


-5.0 


SI 


16818 


241.2 


1 


3.4 


X 


10^ 


R 


1 


Ark 120 (Mrk 1095) 


so/a 


0.0 


SI 


9682 


134.6 


1 


1.84X 


10^ 


R 


1 


Arp 102B 


EO 


-5.0 


L1.8 


7245 


99.7 


1 


2.2 


X 


10^ 


G 


2 


Circinus 


Sb: 


3.0 


S2 


449 


4.0 


2 


1.3 


X 


10« 


M 


3 


Fairall 9 


S? 




SI 


14095 


199.8 


1 


8.0 


X 


10^ 


R 


1 


IC 342 


SABcd 


6.0 


H 


31 


1.8 


3 


<5.0 


X 


10^5 


S 


4 


IC 1459 


E3 


-5.0 


L2 


1691 


29.2 


4 


3.7 


X 


10** 


G 


5 


IC 4329A 


S0+ 


-1.0 


SI 


4813 


65.5 


1 


5.0 


X 


10« 


R 


1 


Milky Way 


Sbc 


3.0 






0.008 


5 


2.95X 


106 


S 


6 


Mrk 79 (UGC 3973) 


SBb 


3.0 


S1.5 


6652 


91.3 


1 


5.2 


X 


10^ 


R 


1 


Mrk 110 


Pair? 




SI 


10580 


147.7 


1 


5.6 


X 


10« 


R 


1 


Mrk 279 (UGC 8823) 


SO 


-1.0 


S1.5 


9129 


126.6 


1 


4.2 


X 


10^ 


R 


1 


Mrk 335 


so/a 


0.0 


Sl.O 


7730 


106.6 


1 


6.3 


X 


10« 


R 


7 


Mrk 509 


comp 




SI 


10312 


143.8 


1 


5.78X 


10^ 


R 


1 


Mrk 590 (NGC 863) 


Sa: 


1.0 


S1.2 


7910 


109.2 


1 


1.78X 


10^ 


R 


1 


Mrk 817 (UGC 9412) 


S? 




S1.5 


9430 


131.0 


1 


4.4 


X 


10^ 


R 


1 


NGC 205 (MllO) 


dE5 


-5.0 


A 


-241 


0.74 


6 


<9.3 


X 


10^ 


S 


8 


NGC 221 (M32) 


E2 


-6.0 


A 


-145 


0.81 


4 


3.9 


X 


106 


S 


9 


NGC 224 (M31) 


Sb 


3.0 


A 


-300 


0.76 


4 


3.3 


X 


10^ 


s 


10 


NGC 598 (M33) 


Scd 


6.0 


H 


-179 


0.87 


6 


<1.5 


X 


10-5 


s 


11 


NGC 821 


E6? 


-5.0 


A 


1735 


24.1 


4 


5.0 


X 


10^ 


s 


12 


NGC 1023 


SBO- 


-3.0 


A 


637 


11.4 


4 


3.9 


X 


10^ 


s 


13 


NGC 1068 (M77) 


Sb 


3.0 


S1.9 


1137 


14.4 


7 


1.6 


X 


10^ 


M 


14 


NGC 2778 


E 


-5.0 




2049 


22.9 


4 


2.0 


X 


10^ 


s 


12 


NGC 2787 


SBO+ 


-1.0 


LI. 9 


696 


7.5 


4 


3.9 


X 


10^ 


G 


15 


NGC 3031 (M81) 


Sab 


2.0 


S1.5 


-34 


3.9 


4 


6.3 


X 


10^ 


S 


16 


NGC 3115 


SO- 


-3.0 


A 


720 


9.7 


4 


9.1 


X 


10« 


s 


17 


NGC 3227 


SABa 


1.0 


S1.5 


1157 


20.6 


7 


3.9 


X 


10^ 


R 


1 


NGC 3245 


SO? 


-2.0 


T2 


1358 


20.9 


4 


2.1 


X 


10^ 


G 


18 


NGC 3377 


E5+ 


-5.0 


A 


665 


11.2 


4 


1.0 


X 


108 


S 


12 


NGC 3379 (M105) 


El 


-5.0 


L2/T2: 


911 


10.6 


4 


1.0 


X 


10« 


s 


19 


NGC 3384 


SBO-: 


-3.0 


A 


704 


11.6 


4 


1.8 


X 


10^ 


s 


12 


NGC 3516 


SBO: 


-2.0 


S1.2 


2649 


38.9 


7 


2.3 


X 


10^ 


R 


7 


NGC 3608 


E2 


-5.0 


L2/S2: 


1253 


22.9 


4 


1.1 


X 


10^ 


s 


12 


NGC 3783 


SBa 


1.0 


SI 


2917 


38.5 


7 


9.4 


X 


10« 


R 


1 


NGC 3998 


SO? 


-2.0 


LI. 9 


1040 


14.1 


4 


5.6 


X 


10« 


S 


16 


NGC 4051 


SABbc 


4.0 


S1.2 


725 


17.0 


7 


1.3 


X 


106 


R 


1 


NGC 4151 


SABab 


2.0 


S1.5 


995 


20.3 


7 


1.53X 


10^ 


R 


1 


NGC 4203 


SABO 


-3.0 


LI. 9 


1086 


14.1 


4 


<1.2 


X 


to^ 


G 


15 


NGC 4258 (M106) 


SABbc 


4.0 


S1.9 


448 


7.3 


4 


4.1 


X 


lO'^ 


M 


20 


NGC 4261 (3C 270) 


E2+ 


-5.0 


L2 


2238 


31.6 


4 


5.2 


X 


10« 


G 


21 


NGC 4291 


E 


-5.0 


A 


1757 


26.2 


4 


1.5 


X 


10^ 


S 


12 


NGC 4342 


SO- 


-3.0 




751 


16.8 


7 


3.4 


X 


108 


S 


22 


NGC 4374 (M84, 3C 272.1) 


El 


-5.0 


L2 


1060 


18.4 


4 


1.6 


X 


109 


G 


23 


NGC 4395 


Sm: 


9.0 


S1.5 


319 


3.6 


7 


<1.1 


X 


10^ 


S 


24 


NGC 4459 


S0+ 


-1.0 


T2: 


1210 


16.1 


4 


6.5 


X 


10^ 


G 


15 


NGC 4473 


E5 


-5.0 


A 


2244 


15.7 


4 


1.0 


X 


108 


S 


12 


NGC 4486 (M87, 3C 274) 


E0+ 


-4.0 


L2 


1307 


16.1 


4 


3.4 


X 


10" 


G 
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TABLE 1: Galaxies with Black Hole Masses — Continued 



3 



Galaxy 


Hubble 


T 


Spectral 


cz 


D 


Ref. 




Method 


Ref. 


Name 


Type 




Class 


(km s^^) 


(Mpc) 




(Mq) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


NGC 4486B 


EO 


-6.0 




1555 


16.1 


8 


6.0 xlO*^ 


S 


26 


NGC 4564 


E 


-5.0 


A 


1142 


15.0 


4 


5.7 xlO^ 


S 


12 


NGC 4593 


SBb 


3.0 


SI 


2698 


39.5 


7 


8.1 xlO*' 


R 


7 


NGC 4594 (M104) 


Sa 


1.0 


L2 


1024 


9.8 


4 


1.1 xlO" 


S 


27 


NGC 4596 


SBO 


-1.0 


L2:: 


1874 


16.8 


7 


5.8 xlO'^ 


G 


15 


NGC 4649 (M60) 


E2 


-5.0 


A 


1117 


16.8 


4 


2.0 xlO" 


S 


12 


NGC 4697 


E6 


-5.0 




1241 


11.7 


4 


1.2 xlO** 


s 


12 


NGC 4945 


SBcd 


6.0 


S2 


560 


4.2 


9 


1.1 xlO*' 


M 


28 


NGC 5548 


so/a 


0.0 


S1.5 


5149 


70.2 


1 


1.23x10^ 


R 


1 


NGC 5845 


E 


-5.0 




1456 


25.9 


4 


3.2 xlO** 


S 


12 


NGC 6251 


EO 


-5.0 


S2 


6900 


94.8 


1 


5.4 xlO^ 


G 


29 


NGC 7052 


E4 


-5.0 




4672 


63.6 


1 


3.6 xlO^ 


G 


30 


NGC 7457 


SO-? 


-3.0 


A 


812 


13.2 




3.4 xlO^ 


S 


12 


NGC 7469 


SABa 


1.0 


Sl.O 


4892 


66.6 


1 


6.5 xlO'' 


R 


1 


PG 0026+129 






QSO 


0.142 


627.4 


1 


5.4 xlO^ 


R 


1 


PG 0052+251 


S 




QSO 


0.155 


690.4 


1 


2.2 xlO^ 


R 


1 


PG 0804+761 






QSO 


0.100 


429.9 


1 


1.89x10^ 


R 


1 


PG 0844+349 






QSO 


0.064 


268.4 


1 


2.16x10^^ 


R 


1 


PG 0953+414 


S 




QSO 


0.239 


1118 


1 


1.84x10* 


R 


1 


PG 1211+143 






QSO 


0.085 


361.7 


1 


4.05x10'' 


R 


1 


PG 1226+023 (3C 273) 


E 




QSO 


0.158 


705.1 


) 


5.5 xlO** 


R 


) 


PG 1229+204 


S 




QSO 


0.064 


268.4 




7.5 xlO^ 


R 




PG 1307+085 


E 




QSO 


0.155 


690.4 




2.8 xlO^ 


R 




PG 1351+640 






QSO 


0.087 


370.7 




4.6 xlO^ 


R 




PG 1411+442 






QSO 


0.089 


379.8 




8.0 xlO'^ 


R 




PG 1426+015 (Mrk 1383) 






QSO 


0.086 


366.2 




4.7 xlO^ 


R 




PG 1613+658 (Mrk 876) 






QSO 


0.129 


565.3 




2.41x10^ 


R 




PG 1617+175 (Mrk 877) 






QSO 


0.114 


494.7 




2.73x10** 


R 




PG 1700+518 






QSO 


0.292 


1406 




6 xlO^ 


R 




PG 1704+608 (3C 351) 


E 




QSO 


0.371 


1857 




3.7 xlO'^ 


R 




PG 2130+099 






QSO 


0.061 


255.3 




1.44x10^ 


R 





Note. — Col. (1) Galaxy name. Col. (2) Revised Hubble type from de Vaucouleurs et al. 1991 (RC3), except for QSOs, 
which is estimated by Hamilton, Casertano, & Turnshek 2001. Col. (3) Morphological type index from the RC3. Col. (4) 
Spectral class of the nucleus from Ho et al. 1997, and otherwise from Whittle 1992a and NED, where A = absorption-line 
nucleus, H = H ii nucleus, L = LINER, S = Seyfert, T = "transition object" (LINER/H ii), 1 = type 1, 2 = type 2. and a 
fractional number between 1 and 2 denotes various intermediate types; uncertain and highly uncertain classifications are 
followed by a single and double colon, respectively. Col. (5) Heliocentric radial velocity (redshift for QSOs) from NED. 
Col. (6) Adopted distance. Col. (7) Reference for D. Col. (8) Black hole mass, scaled to our adopted distances. Col. (9) 
Method for determining Mbu- G, gas kinematics; M, maser kinematics; R, reverberation mapping; S, stellar kinematics. 
Col. (10) Reference for Mbh- 

References. — Distance: (1) Luminosity distance derived from heliocentric redshift, Hq = 75 km s^^ Mpc^^, Om = 
0.3, and Qx = 0.7; (2) Freeman ct al. 1977; (3) McCall 1989; (4) Tonry et al. 2001; (5) Reid 1993; (6) Ferrarese et al. 
2000; (7) TuUy 1988, who also uses our value of Ho; (8) Assumed to be at the distance of NGC 4486; (9) Assumed to be 
at the distance of NGC 5128, which is known from Tonry et al. 2001. 

References. — Black hole mass: (1) Kaspi et al. 2000; (2) Newman et al. 1997; (3) Greenhill et al. 2000; (4) Boker, 
van der Marel, & Vacca 1999; (5) Verdoes Kleijn et al. 2000; (6) Genzel et al. 2000; (7) Ho 1999a; (8) .Jones et al. 1996; 
(9) van der Marel et al. 1998; (10) Kormendy & Bender 1999; (11) Gebhardt et al. 2001; (12) Gebhardt et al. 2000a; (13) 
Bower et al. 2001a; (14) Greenhill et al. 1996; (15) Sarzi et al. 2001; (16) Bower et al. 2001b; (17) Emsellem, Dejonghe, 
& Bacon 1999; (18) Barth et al. 2001; (19) Gebhardt et al. 2000b; (20) Miyoshi et al. 1995; (21) Ferrarese, Ford, & Jaffe 
1996; (22) Cretton & van den Bosch 1999; (23) Bower et al. 1998; (24) Filippenko & Ho 2001; (25) Macchetto et al. 1997; 
(26) Kormendy et al. 1997a; (27) Kormendy et al. 1997b; (28) Greenhill, Moran, & Herrnstein 1997; (29) Ferrarese & 
Ford 1999; (30) van der Marel & van den Bosch 1998. 
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0.13 
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NGC 4945 
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-19.58 


21.81 
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21.59 


25 


39.97 


13 


-14.40 


2.09 



RADIO EMISSION AND BLACK HOLE MASS 5 



TABLE 2: Photometric Parameters — Continued 



Galaxy Name 


Ab 


If" 


log Pe.tot 


Ref. 


log Pe.nuc 


Ref. 


log L-HB 


Ref. 




logfl' 




(mag) 


(mag) 


(W Hz"^) 




(W Hz-1) 




(erg s"^) 




(mag) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


NGC 5548 


0.088 


-21.02 


21.91 


8 


21.86 


17 


41.72 


3 


-19.58 


0.29 


NGC 5845 


0.23 


-18.80 


<19.76 


3 


<19.60 


15 




10 






NGC 6251 


0.38 


-21.40 


24.18 


2 


24.22 


26 


38.79 


14 


-10.90 


6.13 


NGC 7052 


0.52 


-20.73 


22.77 


1 


22.64 


27 


39.22 


15 


-12.18 


4.03 


NGC 7457 


0.22 


-18.73 


<19.18 


3 






<36.82 


6 


<-5.07 




NGC 7469 


0.29 


-21.41 


22.58 


8 


22.47 


17 


41.81 


3 


-19.84 


0.79 


PG 0026+129 


0.31 


-24.35 


23.38 


11 


23.38 




43.42 


16 


-24.61 


-0.19 


PG 0052+251 


0.21 


-23.99 


22.63 


11 


22.63 




43.40 


16 


-24.55 


-0.93 


PG 0804+761 


0.15 


-23.17 


22.72 


11 


22.72 


1 


43.21 


16 


-23.99 


-0.61 


PG 0844+349 


0.16 


-23.30 


21.43 


11 


21.43 




43.07 


16 


-23.57 


-1.74 


PG 0953+414 


0.054 


-25.24 


23.45 


11 


23.45 




44.17 


16 


-26.83 


-1.01 


PG 1211+143 


0.15 


-23.31 


24.39 


11 


24.39 




43.11 


16 


-23.09 


1.18 


PG 1226+023 


0.089 


-20.47 


27.34 


11 


27.34 




44.51 


16 


-27.83 


2.48 


PG 1229+204 


0.12 


-22.61 


21.76 


11 


21.76 




42.95 


16 


-23.22 


-1.26 


PG 1307+085 


0.15 


-24.07 


22.30 


11 


22.30 




43.61 


16 


-25.17 


-1.51 


PG 1351+640 


0.088 


-22.52 


23.34 


11 


23.34 




42.65 


16 


-22.33 


0.67 


PG 1411+442 


0.036 


-22.95 


22.02 


11 


22.02 




43.08 


16 


-23.60 


-1.15 


PG 1426+015 


0.14 


-22.91 


22.29 


11 


22.29 




42.98 


16 


-23.31 


-0.77 


PG 1613+658 


0.11 


-23.43 


23.06 


11 


23.06 




43.29 


16 


-24.22 


-0.36 


PG 1617+175 


0.18 


-23.06 


22.50 


11 


22.50 




43.16 


16 


-23.84 


-0.76 


PG 1700+518 


0.15 


-25.46 


24.23 


11 


24.23 




43.79 


16 


-25.70 


0.22 


PG 1704+608 


0.097 


-25.50 


26.71 


11 


26.71 




43.52 


16 


-24.90 


3.01 


PG 2130+099 


0.19 


-22.61 


22.20 


11 


22.20 




42.95 


16 


-23.22 


-0.82 



Note. — Col. (1) Galaxy name. Col. (2) Galactic extinction in the B band, based on the DTRBE/ IRAS maps of Schlegel, 
Finkbeiner, &: Davis 1998 and the extinction law of Cardelli, Clayton, &: Mathis 1989. Col. (3) Total absolute B magnitude, 
corrected for Galactic extinction; data for nearby galaxies from RC3 or LEDA; data for PG QSOs from Schmidt &: Green 1983; 
data for Milky Way assuming Mk = -24.06 mag (Malhotra et al. 1996) and B - K = 3.53 mag (for T = 3; de Jong 1996). Col. 
(4) Logarithm of the total 6 cm spectral power obtained from low-resolution (beam >l' in most cases) observations; in a few cases 
where 6 cm data are unavailable we extrapolated measurements from other wavelengths assuming oc v^"'" , as found by Ho & 
Ulvestad 2001. Col. (5) Reference for Pe.tot- Col. (0) Logarithm of the nuclear 6 cm spectral power. Col. (7) Reference for Pe.nuc. 
Col. (8) Logarithm of the Hd (narrow + broad components) luminosity, corrected for Galactic extinction. Col. (9) Reference 
for Lh,3- Col. (10) Absolute B magnitude of the nucleus, obtained from the Lh/3 — Mb relation of Ho & Peng 2001. Col. (11) 
Logarithm of the ratio of the nuclear luminosities in the radio to the optical B band, defined as R' = L^(6 cm) /L'^{B), with L'^{B) 
computed from M'g . 

References. — Integrated radio emission: (1) Becker, White, &: Edwards 1991, resolution 3'.5; (2) White &: Becker 1992, resolution 
12', extrapolated from 20 cm; (3) Condon et al. 1998, resolution 0'.75, extrapolated from 20 cm; (4) Gregory et al. 1994, resolution 
4'.2; (5) Whittle 1992a, extrapolated from 20 cm; (6) Wright et al. 1996, resolution 4'.2; (7) Berkhuijsen 1984, extrapolated from 408 
MHz assuming /„ oc f"" '^; (8) Rush, Malkan, & Edelson 1996, resolution 1'.5; (9) Condon 1987, resolution 1', extrapolated from 
20 cm; (10) Griffith et al. 1994, resolution 4'.2; (11) Kellermann et al. 1989. 

References. — Nuclear radio emission: (1) Assumed to be equal to the total radio power; (2) Barvainis, Lonsdale, &: Antonucci 
1996; (3) Puschell et al. 1986; (4) Elmouttie et al. 1998; (5) Turner & Ho 1983; (6) Sadler, Jenkins, & Kotanyi 1989; (7) Unger et 
al. 1987, extrapolated from 20 cm; (8) Ekers et al. 1983, for Sgr A* only; (9) Ulvestad & Wilson 1984a; (10) Kukula et al. 1995, 
extrapolated from 3.6 cm; (11) Singh & Westergaard 1992; (12) Heckman, Balick, & Crane 1980; (13) Crane, Dickel, & Cowan 
1992, extrapolated from 3.6 cm; (14) Hummel et al. 1987, extrapolated from 20 cm; (15) Wrobel & Heeschen 1991; (16) Ulvestad 
fc Wilson 1984b; (17) Ho & Ulvestad 2001; (18) Fabbiano, Gioia, & Trinchieri 1989; (19) Wrobel & Heeschen 1984; (20) Jenkins, 
Pooley, &: Riley 1977; (21) Biretta, Stern, & Harris 1991; (22) Hummel, van der Hulst, &: Dickey 1984; (23) Spencer &: Junor 1986; 
(24) Birkinshaw &: Davies 1985; (25) Elmouttie et al. 1997; (26) Jones et al. 1986; (27) Morganti et al. 1987. 

References.— Nuclear H/? emission: (1) Tadhunter et al. 1993; (2) Oke k Goodrich 1981; (3) Whittle 1992a; (4) Halpern et 
al. 1996; (5) Oliva et al. 1994, corrected for extinction; (6) Ho, Filippenko, & Sargent 1997; (7) Verdoes Kleijn et al. 2000; (8) 
extrapolated from the ix(2 — lOkeV) — Lua relation of Ho et al. 2001, based on the Chandra 2 — 10 keV detection of Sgr A* 
reported by Baganoff et al. 2001; (9) Heckman 1996, spatially integrated over the central few hundred pes; (10) no H/9 measurement 
available; optical continuum upper limit (given in col. 11) obtained from archival if^T image, following the method of Ho & Peng 
2001 and Ravindranath et al. 2001; (11) Ho, Filippenko, & Sargent 1990; (12) unpublished, but based on data from Ho et al. 1995, 
1997; (13) Moorwood & Oliva 1988, based on Br7 flux and an intrinsic ratio H/9/Br7 = 36 appropriate for Case B recombination in 
low-density gas at % = lO" K (Osterbrock 1989); (14) Shuder & Osterbrock 1981; (15) Morganti, Ulrich, & Tadhunter 1992 give 
the flux of Ha+[N ll], which we use, along with the [N ii]/Hq; ratio of van den Bosch & van der Marel 1995, to estimate the H/3 
flux; (16) Boroson & Green 1992, obtained from H/3 equivalent widths and optical continuum flux densities (see Ho & Peng 2001). 
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2.1. Black Hole Masses 

Our primary sample of BH masses comes from spatially re- 
solved observations of gas and stellar kinematics. Reviews of 
these methods can be found in Kormendy & Richstone (1995) 
and Ho (1999a). Table 1 represents an update of the compi- 
lation of BH masses given in Ho (1999a), supplemented with 
newly published values. The majority of the gas measure- 
ments consist of optical observations of ionized gaseous disks, 
all done with the Hubble Space Telescope (HST); four galaxies 
(Circinus, NGC 1068, NGC 4258, and NGC 4945) exploit the 
availability of strong water maser emission to probe the nuclear 
kinematics using radio interferometry. The broad-line radio 
galaxy Arp 102B presents a special case in which a periodic sig- 
nature in its emission-Une light curve has been interpreted as the 
orbital period of an accretion disk (Newman et al. 1997). With 
the exception of Sgr A* in the Galactic Center, whose mass has 
been determined through proper motions and radial velocities 
of individual stars (Eckart & Genzel 1997; Ghez et al. 1998; 
Genzel et al. 2000), all the stellar-based masses come from in- 
tegrated spectroscopy, a few from the ground, but largely from 
HST. 

Magorrian et al. (1998) published Mbh for a significant num- 
ber of nearby early-type galaxies based on axisymmetric, two- 
integral dynamical modeling of HST images and ground-based 
stellar spectroscopy. We do not use these masses, however, be- 
cause they are subject to considerable uncertainties due to the 
simplified modehng and the low spatial resolution of the spec- 
tra. As commented by a number of authors (van der Marel 
1999; Ferrarese & Merritt 2000; Gebhardt et al. 2000b), the 
two-integral assumption will tend to overestimate Mbh- We 
illustrate this explicitly for the 16 galaxies in Magorrian et 
al.'s sample that overlap with the sample assembled in Table 1 
(Fig. 1). The average offset for these objects is a factor of 3.3. 
Figure lb indicates that the discrepancy seems to be slightly 
more severe for "core" galaxies than for "power-law" galaxies, 
systems with and without a resolved break in their inner surface 
brightness profiles, respectively (Lauer et al. 1995; Faber et al. 
1997). As a class, core galaxies tend to be luminous, pressure- 
supported systems with boxy isophotes (e.g., Faber et al. 1997; 
Ravindranath et al. 2001), precisely those most prone to having 
an anisotropic velocity distribution and thus most ill-suited for 
two-integral modeling. Figure lb offers mild support for this 
picture. 

We include several meaningful upper limits. Four come from 
stellar-kinematical constraints placed on the central nuclear star 
cluster (IC 342, M33, NGC 205, and NGC 4395), and a fifth 
derives from HST gas kinematics (NGC 4203). 

Direct dynamical measurements become unfeasible for lu- 
minous or distant AGNs. The bright continuum emission of the 
active nucleus nearly always overpowers the stellar absorption 
lines near the center, and in many cases the extended line emis- 
sion can be strongly perturbed by nongravitational forces. For 
these objects we must rely on more indirect methods to estimate 
the central masses. A promising technique employs "reverber- 
ation mapping" (Blandford & McKee 1982) to determine the 
size of the broad-line region, r, which when combined with a 
characteristic velocity dispersion of the line-emitting gas, v, as 
reflected in the observed line widths, yields an estimate of the 

virial mass, Mvir = rv^/G. Choosing v = ^FWHM(H/3) for 
random, isotropic orbits, several studies have computed virial 
masses in this fashion for Seyfert 1 nuclei (Ho 1999a; Wandel, 
Peterson, & Malkan 1999) and low-redshift quasars (Kaspi et 



al. 2000). Ho (1999a) used FWHM(H/3) from single-epoch 
spectra, whereas Wandel et al. (1999) argue that the variable 
component (rms) of the spectrum should yield a more faith- 
ful representation of the velocity field associated with the time 
lag used to calculate r. As Wandel et al. note, however, the 
simple expectation that time-averaged spectra should yield nar- 
rower line widths than rms spectra appear not to hold in all 
objects. In practice, it seems difficult to justify one choice over 
the other (Kaspi et al. 2000). A more serious uncertainty lies 
in the choice of v. McLure & Dunlop (2001) advocate that 
a disk component with v = |FWHM(H/?) yields Mbh values 
for AGNs that agree better with the MBH-^-buige relation for 
weakly active galaxies, thereby obviating the apparent discrep- 
ancy found by Ho (1999a) and Wandel (1999). The situation 
is far from clear, however, since AGN masses computed as- 
suming random orbits agree, to first order, with the MBH-Ce 
relation (Gebhardt et al. 2000c; Nelson 2000; Ferrarese et al. 
2001). Notwithstanding this and other possible systematic un- 
certainties (Krolik 2001), reverberation mapping appears capa- 
ble of delivering Mbh for AGNs with an accuracy of a factor of 
~2-3. 

Table 1 lists Mbh based on reverberation mapping for 20 
Seyfert 1 galaxies and 17 low-redshift quasars. For the sake 
of uniformity, we use the compilation of Kaspi et al. (2000) 
for the majority of the objects; a few not included in that study 
were taken from Ho (1999a). 

Galaxy distances enter linearly into the mass determination 
for all methods except reverberation mapping, wherein the size 
scale depends only on the hght-travel time between the central 
continuum source and the line-emitting gas. Thus, it is impor- 
tant to pay close attention to the reUabiUty of the adopted dis- 
tances. Whenever possible we use Tonry et al.'s (2001) homo- 
geneous database of distances based on /-band surface bright- 
ness fluctuations. 

2.2. Nuclear Luminosities 

Our subsequent analysis primarily will examine the connec- 
tion between Mbh and two measures of the radio output, namely 
the absolute spectral power and the radio-to-optical luminosity 
ratio. Since we are interested in quantities pertaining to the nu- 
cleus, we take great effort to assemble nuclear radio and optical 
data. As emphasized in several recent studies of low-luminosity 
galactic nuclei (Ho 1999c; Ho et al. 2000, 2001; Ho & Peng 
2001), high angular resolution is of paramount importance for 
isolating the central emission from the surrounding galaxy. 

For the radio band, we make use of interferometric data ob- 
tained at 6 cm with beam sizes ^5". The nuclear 6 cm spectral 
power, /'6,nuc, computed from the observed flux densities as- 
suming isotropic emission, represents the integrated emission 
from all components considered associated with the "active" 
(nonstellar) nucleus. In the case of the brighter AGNs, this of- 
ten includes some extended structures, such as jet- like linear 
features, in addition to the central core. For comparison with 
recent results from the literature, we have also collected radio 
data for the integrated emission from the whole galaxy (host 
plus nucleus), /^6,tot> which we approximate with low-resolution 
(beam ;^1') measurements. A minority of the data were ac- 
quired at wavelengths other than 6 cm, and these were extrap- 
olated to 6 cm assuming oc u~^-^ , the median spectrum be- 
tween 6 and 20 cm for Seyfert nuclei found by Ho & Ulvestad 
(2001). 

Ho & Peng (2001) demonstrated the utility of //5r images 



RADIO EMISSION AND BLACK HOLE MASS 



7 



10 



9 - 



"T — I — I — I — I — I — I — I — I — I — r 



"T — I — I — I — I — r 



• stellar kin. 
i gas kin. 



I I I I I 



(a) 



I I I I I I I 



7 8 9 

log [MBH(Others) / M^] 



10 



0.5 



o 







-0.5 



• * core profile 

OA power-law profile 



I 



• o 




-16 -18 -20 -22 

M° (mag) 



Fig. 1. — (a) Comparison between black hole masses from Magorrian et al. (1998) with those obtained from others. The masses from Magorrian et al. are 
systematically higher (h) The difference in log Mbh as a function of the luminosity and surface brightness profile of the galaxy. The solid line denotes equality. 
The most discrepant masses are found in luminous galaxies with core-type profiles. 



for separating the optical continuum cores in a sample of nearby 
Seyfert 1 galaxies. Although //primages are available for most 
of the weakly active galaxies, here we take a different approach. 
A number of the objects have type 2 nuclei\ which, according 
to the AGN unification picture (Antonucci 1993), implies that 
the nuclear continuum should be hidden from direct view, or at 
least appreciably extincted. Instead, we estimate the continuum 
strength indirectly through the known (Yee 1980; Shuder 1981) 
correlation between H/3 luminosity and B-band absolute mag- 
nitude for type 1 AGNs, as recently calibrated by Ho & Peng 
(2001). The line luminosity (broad and narrow components 
combined) is, in principle, a more isotropic quantity than the 
optical continuum luminosity. We collected nuclear Balmer- 
line luminosities (or upper limits thereof for sources lacking 
line emission), translated them to H/3 if necessary^ (only for a 
few cases; see Table 2), and then applied the Lh/3 -Mb conver- 
sion (Ho & Peng 2001) to arrive at M'g. We use the notation 
M'g to distinguish it from the directly measured quantity Mb. 
Note that Mb is available for all the quasars (Schmidt & Green 
1983) and for a number of the Seyfert 1 nuclei (Ho & Peng 
2001); for the bright, variable Seyfert nuclei studied with rever- 
beration mapping, it can be ascertained quite reUably from the 
published spectrophotometry. For internal consistency, how- 
ever, we follow the same procedure as adopted for the weakly 
active sample. 

As discussed in Ho & Peng (2001), the individual Une lumi- 
nosities can be quite uncertain, and the low-luminosity end of 
the Lup -Mb relation has significant scatter. We do not expect 
the continuum magnitudes to be very accurate for any given 
object, but it is hoped that the statistical results are more robust. 

We use the spectral radio luminosity and the optical nuclear 
luminosity to calculate the equivalent of the standard radio- 

' Type 1 and type 2 AGNs are defined as those with and without detectable broad 
^ We adopt an intrinsic ratio of Ha/H^ = 3.1, as might be appropriate for the phys 



loudness parameter/?, R' = L^(6 cm) /Ll(B). Following com- 
mon practice (Visnovsky et al. 1992; Stocke et al. 1992; Keller- 
mann et al. 1994), we set the boundary between "radio-loud" 
and "radio-quiet" classes at R' = 10. 

2.3. The Special Case of the Galactic Center 

The Galactic Center warrants some individual attention. Be- 
cause of its proximity, the "nucleus" of the Milky Way — iden- 
tified with Sgr A* — is overresolved with respect to other 
galaxies, thus rendering comparisons somewhat ambiguous. 
In our subsequent discussions, we relax the definition of the 
Galactic Center to include increasingly larger areas surround- 
ing Sgr A*. At radio wavelengths, two distinct regions can be 
identified in the vicinity of Sgr A*, namely Sgr A West and 
Sgr A East (Ekers et al. 1983). Each of the latter two compo- 
nents is ~40 times brighter than Sgr A* at 6 cm. To achieve a 
linear scale roughly equivalent to that sampled in most of the 
external galaxies, however, one must extend to dimensions of 
^l°x 1°, or approximately 150 pcxlSO pc. The 6 cm flux 
density on this scale (600 Jy; Mezger & Pauls 1979) is roughly 
10 times higher than the entire Sgr A complex combined. 

Evaluating the radio-loudness parameter for the Galactic 
Center requires knowledge of its intrinsic optical continuum 
luminosity, ideally measured on the various scales discussed 
above for the radio emission. The optical emission, of course, 
is not directly observable because of the tremendous opacity 
along our line of sight. Instead, we proceed as follows. We 
utilize the luminosity measured in the hard X-ray (2-10 keV) 
band, corrected for photoelectric absorption, to estimate the in- 
trinsic (unextincted) Ha luminosity using the relation between 
Lx(2 - 10 keV) andLna proposed by Ho et al. (2001) for nearby 
galactic nuclei. Next, Lh/3 follows straightforwardly from Lna 

ssion lines, respectively. 

conditions in active nuclei (e.g., Gaskell & Ferland 1984). 
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log (M3„ / M3) log (M3H / M3) 

Fig. 2. — Dependence of (a) total and (b) nuclear 6 cm spectral power on black hole mass. The symbols for the different methods of mass determination are given 
in the legend, below which representative error bars are plotted. Arrows indicate upper limits. The clashed line in each graph gives the regression fit proposed by 
Franceschini et al. (1998); the solid lines are discussed in the text. In panel (b) we have labeled the outlier M31. Four values of Pe.nuc are given for the Milky Way 
(MW), connected by the dotted line; in increasing value, they are for Sgr A*, Sgr A West, Sgr A West+East, and the central 1° x 1° (~150pcxl50 pc) (see § 2.3). 



for an assumed Ha/H/3 ratio. And finally, the relation between 
Lh/3 and Mb (Ho & Peng 2001) yields the B-band continuum 
flux density needed to calculate the R parameter 

Baganoff et al. (2001) recently detected Sgr A* unam- 
biguously in the hard X-rays with Chandra; the unabsorbed 
2-10 keV luminosity has a surprisingly low value of 2.4x 10^^ 
erg s"' . The hard X-ray emission on larger scales comes from 
observations performed using ASCA by Koyama et al. (1996). 
On dimensions which encompass Sgr A West, Lx(2- 10 keV) 
= 7.7x10^^ erg s"', while for flie l°x 1° region they find 
Lx(2- 10 keV) = 7.7x10^'' erg s"^ We could not locate a 
definitive value for the X-ray luminosity of Sgr A East. Follow- 
ing the above procedure, we find that Sgr A*, Sgr A West, and 
the l°x 1° region have \ogR = 4.9, 3.8, and 4.0, respectively. 
The Galactic Center, irrespective of one's exact definition of 
its boundaries, evidently is extremely "radio loud" according to 
the conventional R parameterization. 

Admittedly, the above conclusions for the Galactic Center 
depend strongly on the applicability of the conversion factors 
used to translate Lx(2- 10 keV) to Mg, ones which were origi- 
nally derived for more luminous, more active type 1 nuclei (see 
Ho et al. 2001; Ho & Peng 2001). The spectral classification of 
the Galactic Center is unknown. Ho et al. (2001) find that type 2 
nuclei with detectable X-ray cores generally have Lx/Lua val- 
ues that are a factor of ^10 lower than in type 1 objects. Even 
a factor of 10 error, or greater, however, cannot erase the large 
R values given above. And although a direct relation between 
Lh/3 and Mb does not exist for type 2 objects (the optical contin- 
uum is weak or obscured), to the extent that the line emission is 
powered by photoionization and Mb traces the low-energy tail 
of the ionizing continuum, the two quantities should roughly 
scale with one another in type 2 objects as they do in type 1 
systems. 



2.4. Error Estimates 

Systematic errors affect many of the quantities used in this 
paper. Although formal uncertainties are not specified explic- 
itly for the data presented in Tables 1 and 2, we wish to alert 
the reader of their hkely magnitude and impact on our analysis. 
The representative error bars shown in the figures of this paper 
are meant to capture the assessment given in this section. 

Ho & Peng (2001; see § 3.3) give a fairly thorough account 
of the error budget associated with the radio and optical lumi- 
nosities. We will not repeat the details here, except to reiterate 
that the typical uncertainties for the radio powers, optical line 
luminosities, and R' are ^0.2, 0.3, and 0.5 dex, respectively. 

At the moment, the uncertainties on Mbh are still quite var- 
ied. The masses for the BHs in the Galactic Center (Genzel 
et al. 2000) and NGC 4258 (Miyoshi et al. 1995) are known 
with high confidence, on the order of 10%. Masses based on 
three-integral modeling of stellar kinematics are accurate to 
^0.3 dex on average (Gebhardt et al. 2000b), but the few that 
are still based on two-integral models may undergo more sig- 
nificant revisions in the future. While it is generaUy thought 
that gas-kinematical methods are less prone to modeling un- 
certainties than those based on stellar kinematics, under some 
circumstances our inabiUty to treat realistically certain effects 
such as asymmetric drift may lead to serious systematic errors. 
In the case of IC 1459, the effect of asymmetric drift on Mbh 
may be as large as a factor of 4 (Verdoes Kleijn et al. 2000). 
This correction, on the other hand, is much less significant in 
NGC 3245, whose Mbh has been determined to an accuracy of 
~25% (Earth et al. 2001). Lastly, we concluded in § 3.1 that 
reverberation-mapping masses are probably accurate to a factor 
of 2- 3. Thus, an overall uncertainty of 0.3 dex for Mbh seems 
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Fig. 3. — Dependence of total 6 cm spectral power on the absolute B mag- 
nitude of the galaxy. Symbols as in Figure 2. Representative error bars are 
plotted in the upper left comer; arrows indicate upper limits. The solid line 
shows the formal linear regression fit. 

reasonable. 

3. PARAMETER CORRELATIONS 
3.1. Radio Power and Black Hole Mass 

Concentrating first on the purported dependence of total ra- 
dio power on Mbh, Figure 2a shows that these two quantities 
are not well correlated for the sample of weakly active galax- 
ies, contrary to the results of Franceschini et al. (1998). Laor 
(2000) recently arrived at a similar conclusion using the less 
rehable (see § 2.1) BH masses from Magorrian et al. (1998). 
Although the best-fit regression Une proposed by Franceschini 
et al. (Lrad oc Mg^; dashed line) does roughly follow the trend 
of the data, the scatter about the fit is enormous. The distribu- 
tion of points may be tracing an upper envelope. The absence 
of points on the upper left corner of the diagram is real: lumi- 
nous radio sources associated with low-mass BHs (which pref- 
erentially inhabit low-mass bulges) can hardly be missed. The 
lower right portion of the graph, on the other hand, is mostly 
populated by upper hmits, and the remaining blank region may 
reflect the observational bias against finding faint radio sources 
in the most luminous, on average more remote, galaxies which 
house the heftiest BHs. 

The form of the ridge-line of the upper envelope, which ap- 
proximately follows Lrad OC ^BH~^'^> itself cau be explained as 
an indirect by-product of two known, more fundamental corre- 
lations. Although the physical cause is not well understood, it 
has long been known that the integrated radio emission of early- 
type galaxies increases with their total optical luminosity (e.g., 
Auriemma et al. 1977; Fabbiano, Gioia, & Trinchieri 1989; 
Sadler, Jenkins, & Kotanyi 1989; Calvani, Fasano, & Frances- 
chini 1989; Ho 1999b). The total mass (Heckman 1983) or 
pressure (Whittle 1992b) of the bulge component has been sug- 
gested as a parameter that could affect the efficacy of generating 
radio emission. In previous studies, the relation can be 



Fig. 4. — Correlation between nuclear 6 cm spectral power and H/3 (broad + 
narrow components) luminosity. Symbols as in Figure 2. Representative error 
bars are plotted in the upper left comer; arrows indicate upper limits. We have 
labeled tiie outlier M31. Four values of Pj.nuc are given for the Milky Way 
(MW), connected by the dotted line; in increasing value, they are for Sgr A*, 
SgrAWest, Sgr A West+East, and the central l°x 1° (~150pcxl50pc)(see 
§ 2.3). 

described as Lr^d oc ifp,, generally with j3 « 1-3. Figure 3 
shows the strong correlation between total radio emission and 
integrated absolute B magnitude for our sample. The gener- 
alized KendaU's r test (Isobe, Feigelson, & Nelson 1986) re- 
turns a correlation coefficient of -1.1 at a significance level 
>99.99%. A linear regression fit {solid line) using Schmitt's 
(1985) method gives logPe.tot oc -0.99M2^, or Lrad oc Ll^^. 
Now, since Mbh scales roughly linearly with the optical lu- 
minosity of the bulge (see references in § 1), which in bulge- 
dominated systems is comparable to the integrated light of the 
galaxy, it follows that Lrad oc Mgy, as observed (Fig. 2a; solid 
line). 

The active galaxies in Figure 2a do form a fairly well de- 
fined correlation, but this is simply a manifestation of the fact 
that for these objects the AGN component dominates the inte- 
grated emission, and, as we now argue, the relation between the 
nuclear radio power and Mbh is physically meaningful. 

The distribution of Pe.nuc vs. Mbh also resembles an upper 
envelope (Fig. 2b), whose overall appearance is similar to that 
of the Pg tot vs. Mbh diagram. Again, the weakly active galax- 
ies do not trace the Franceschini et al. relation. The location 
of the points for Sgr A* and the nucleus of M31 are particu- 
larly striking. They are the only objects with stellar-kinematical 
masses in the range Mbh ~ 10^ to few x 10^ Mq that have de- 
tected radio cores. Both are extremely weak, with /^6,nuc ^10'^ 
W Hz~\ approximately 3 orders of magnitude lower than the 
upper limits for the more distant galaxies. The position of the 
Galactic Center point, however, would migrate upwards (see 
points joined by dotted line) depending on one's definition of 
the "center," in which case M31 would be a distinct outUer. In 
any case, we suspect that the undetected galaxies could have 
tiny radio cores like Sgr A* or the nucleus of M3 1 if they were 
to be observed at much higher sensitivity and resolution. Al- 
though the distribution of points in the lower right comer of 
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Fig. 5. — Distribution of tlie nuclear radio-to-optical luminosity ratio R' (see text) versus (a) black hole mass and (b) A = Lboi/^E- The solid line marks the 
formal division between radio-loud and radio-quiet objects, i?' = 10. The dashed line in panel {b) is the best-fitting linear regression line. The legends indicate 
the meaning of the symbols. Representative error bars are plotted in the lower left comer of each panel; arrows indicate limits. The most massive BHs tend to be 
exceptionally radio-loud, highly sub-Eddington systems. We have labeled the outlier M31. Three points are plotted for the Milky Way (MW), coimected by the 
dotted line; in increasing value of R' , they are for Sgr A West, the central 1° x 1° (--jISO pcx 150 pc), and Sgr A* (see § 2.3). 



Figure 2b remains to be determined, it is likely that Mbh will 
continue to be a very poor predictor of radio power. Consider 
M3 1 and M8 1 . Their BH masses differ by only a factor of ~2, 
but the contrast in their radio core powers is a factor of nearly 
10^ 

The ridge-line of the upper envelope, on the other hand, is 
fairly well defined, and as in Figure 2a, the "zone of avoidance" 
in the upper left corner is genuine. Ignoring for the moment the 
handful of high-power sources with nue > 10^ W Hz"\ it 
appears that at a given mass, the maximum core radio power 
increases roughly linearly with Mbh- A natural explanation for 
the functional form of the ridge-line can be found by appeal- 
ing to Eddington-hmited accretion and recognizing that radio 
power also scales roughly linearly with accretion (optical) lumi- 
nosity for relatively active AGNs (Ho & Peng 2001), or equiv- 
alently, with mass accretion rate. It is worth remarking that the 
sample of the present study, although significantly more het- 
erogeneous and diverse than that investigated by Ho & Peng 
(2001), also displays a reasonably strong correlation between 
radio and optical luminosity (Fig. 4)^. The maximum lumi- 
nosity output of an accreting object in an isotropic, homoge- 
neous system is set by the Eddington luminosity, Le = 1.3 x 
Kfi^iM^ii/ Mr.) erg s"'. Setting Lboi = ^e, and approximating 
the bolometric luminosity with Lboi ~ cbLb ~ CBCrPe,nuc, we 
immediately arrive at Pe.nuc oc Mg^. The constant cg ranges 
from ~ 1 1 - 17 for luminous AGNs and quasars (Sanders et al. 
1989; Elvis et al. 1994) to ~24 for very low-luminosity, possi- 
bly advection-dominated systems (median value for the 12 ob- 
jects studied by Ho 1999c and Ho et al. 2000). The constant 
Cr can be obtained from the radio-optical continuum correla- 
tion of Ho & Peng (2001), and it differs slightly for radio-loud 

^ As with the results reported in Ho & Peng (2001), we have confirmed that this 
statistical method used to evaluate partial correlations with a third variable. 



compared to radio-quiet objects. The solid line in Figure 2b 
was obtained by choosing cb = 17 (Sanders et al. 1989) and 
the radio-loud branch of the Pe.nuc ^ relation (Ho & Peng 
2001). The agreement between the Une and the boundary of the 
upper envelope is surprisingly good, both for the slope and the 
intercept. Nearly all of the AGNs fall in a band bracketed by 
Lboi/LE ~ 0.01-1. Notably, many of the moderately active but 
lower luminosity AGNs (nearby Seyfert 2 nuclei and LINERs) 
whose masses have been determined by maser or ionized-gas 
kinematics are also broadly distributed among the more lumi- 
nous AGNs. By contrast, the optically and radio quiescent sys- 
tems, which comprise all the objects with stellar-based masses, 
uniformly occupy the highly sub-Eddington regime of the dia- 
gram. A small cluster of radio-luminous objects (e.g., 3C 120, 
3C 273, 3C 351, 3C 390.3) lie above the Lboi = Le line, but this 
is not unexpected because objects with powerful radio jets are 
known to follow a steeper radio-optical correlation (e.g., Ser- 
jeant et al. 1998; Willott et al. 1999). 

3.2. Radio Loudness and Black Hole Mass 

BH accretion in galactic nuclei invariably generates radio 
emission. The fundamental parameters responsible for the 
tremendous range of the observed strength of the radio output, 
however, are not well estabUshed and have been largely a sub- 
ject of speculation. Recent advances in high-resolution imaging 
of quasars consistently suggest that radio-loud objects reside in 
hosts which lie on the top end of the galaxy luminosity function, 
whereas the hosts of radio-quiet sources generally span a wider 
range of luminosities. The degree of radio loudness, therefore, 
seems to depend on galaxy mass. In view of the link between 
BH mass and bulge mass (see § 1), it is reasonable to deduce 

elation is not a spurious distance effect. See Ho & Peng for a discussion of the 
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that radio loudness would depend on Mbh- Laor (2000) exam- 
ined this issue using a sample of low-redshift (z < 0.5) quasars 
from the Palomar-Green (PG) survey (Schmidt & Green 1983). 
The radio properties of the PG sources are known (Kellermann 
et al. 1989), and approximate virial masses can be estimated 
from the H/3 line widths and optical continuum luminosities 
following empirical calibrations derived from AGN variability 
studies (Laor 1998; Kaspi et al. 2000). Laor (2000) concluded 
that the radio-loudness parameter /? is a strong function of Mbh- 
most radio-loud objects (defined by log/? > 1) have Mbh > W 
Mq, whereas nearly all quasars with Mbh < 3x10^ Mq are ra- 
dio quiet. 

The sample considered in this study, which spans a much 
wider gamut of activity level, yields a more complex picture. 
As shown in Figure 5a, the clean segregation in mass between 
radio-loud and radio-quiet objects suggested by Laor disap- 
pears. The radio-loud plane (log/?' > 1) is richly populated 
with masses ranging from log Mbh ~ 9.5 to 6.0, and possibly 
even lower if we consider the upper limits. Indeed, most or 
all objects with log Mbh ^8.5 are radio loud, but radio-loud 
objects are by no means restricted to the high-mass domain. 

To explore other factors which may govern the distribution 
of points on this plot, we have coded the symbols according 
to the bolometric luminosity of the nucleus normalized to the 
Eddington luminosity, A = Lboi/^^E- As in § 3.1, we estimate 
Lboi crudely from Lb (which is based on Lh/s). An intriguing 
pattern emerges. The majority of sources with A ^ 1 land in 
the radio-quiet regime; those with A <; 10"^ fall exclusively in 
radio-loud territory; and objects with intermediate values of A 
straddle the (somewhat arbitrary) R' boundary. (M3 1 is a per- 
sistent outlier.) The dependence of R' on A is shown explictly 
in Figure 5b; the objects form a striking inverse correlation, al- 
beit with substantial scatter. The majority of the strongly active 
nuclei (Seyfert Is and quasars) are characterized by log A ^ 
-2 and log/?' < 1 . With the exception of the deviant point for 
M31, all of the weakly active objects (those with Mbh based on 
spatially resolved kinematics) are confined to log A < -2 and 
log/?' > 1. More quantitatively, the generalized Kendall's r 
test, with M31 omitted, gives a correlation coefficient of -0.97 
at a significance level >99.99%. A hnear regression fit (dashed 
line) using Schmitt's (1985) method, which treats censoring in 
both variables, gives 

log/?' = -(0.50 ± 0.07) logA + (0.27 ±0.24). 

Since A depends on the mass accretion rate M (see, e.g.. Fig. 7 
of Narayan et al. 1998b), an immediate consequence of the 
/?' - A inverse correlation is that the degree of radio loudness 
depends strongly on M. 

Although the strong inverse correlation between /?' and A 
may superficially resemble a mutual dependence of these vari- 
ables with optical luminosity (/?' oc A oc Lg), we note that 
the slope of the /?' - A relation differs significantly from -1 . We 
do not believe that this is the primary driver of the observed 
correlation. 

4. DISCUSSION AND SUMMARY 

A number of statistical studies of luminous AGNs indi- 
cate that the processes responsible for the generation of radio 
sources depend directly on BH mass (McLure et al. 1999; Nel- 
son 2000; Laor 2000; McLure & Dunlop 2001). Lacy et al. 
(2001) suggest that, in addition to the primary dependence on 
Mbh, the radio luminosity is also a weak function of L/L^. 



Along similar lines, Franceschini et al. (1998) proposed that 
the weakly active nuclei in nearby galaxies also obey a corre- 
lation between radio luminosity and Mbh, roughly of the form 
Z-rad oc Mg^. The correlation evidently holds for the core radio 
emission as well as for the integrated radio emission. Frances- 
chini et al. (1998) and Di Matteo et al. (2001) argue that the 
slope of the Lrad -Mbh relation can be explained qualitatively in 
the context of advection-dominated accretion of hot plasma un- 
dergoing Bondi-type inflow. This interpretation, however, has 
been questioned by Yi & Boughn (1999), who noted that, given 
the observed radio powers and the critical assumption of Bondi 
accretion, the ADAF model predicts far more X-ray radiation 
than is actually detected. 

We have reexamined these issues using a comprehensive 
compilation of up-to-date BH masses and photometric param- 
eters. Our sample includes all the nearby, weakly active or 
inactive galaxies which have reliable BH masses determined 
through spatially resolved kinematics, as well as all the AGNs 
(Seyfert 1 nuclei and quasars) for which virial masses have been 
derived through reverberation mapping. Our main results lead 
to the foUowing conclusions. 

1. There is no simple relation between integrated radio lu- 
minosity and Mbh- The distribution of objects is consistent 
with either an upper envelope or a loose correlation of the 
form Lrad C!C Mg^"^^, similar to that found by Franceschini et 
al. (1998), but we offer a different interpretation. We suggest 
that the integrated Lrad -Mbh relation arises indirectly from two 
known, physically more fundamental correlations, namely that 
between integrated radio luminosity and optical bulge luminos- 
ity (mass) and that between bulge luminosity (mass) and BH 
mass. 

2. There is no simple relation between core radio luminos- 
ity and Mbh- The distribution of objects appears to follow an 
upper envelope defined by Eddington-hmited accretion. At any 
given value of Mbh, the maximum luminosity attained is set by 
Lboi = Le, and as discussed by Ho & Peng (2001), the core ra- 
dio power traces the accretion luminosity. The majority of the 
more active nuclei in our sample appear to be characterized by 
Lboi ~ (0.01- 1)Le- Not surprisingly, the more quiescent ob- 
jects are radiating at only a tiny fraction of the Eddington rate. 

3. There is no simple relation between the radio-loudness pa- 
rameter/? and Mbh- Specifically, we do not find the clean divi- 
sion between R and Mbh suggested by Laor (2000). Radio-loud 
nuclei are not confined solely to galaxies with the most massive 
BHs, but instead can inhabit galaxies with Mbh as low as 10^ 
Mq, or perhaps even less. The nucleus of the Milky Way pro- 
vides an interesting, if unfamiliar, illustration. As a corollary 
of this result, radio-loud nuclei are not restricted to early-type 
(SO and elliptical) galaxies, as conventionally thought; rather, 
they can be hosted by galaxies of a variety of morphological 
types, including disk-dominated spirals. Both of these results 
have been foreshadowed by recent investigations of the nuclear 
spectral energy distributions of nearby low-luminosity AGNs 
(Ho 1999c; Ho et al. 2000). Ho & Peng (2001) specifically 
challenged the traditional notion that Seyfert nuclei, the major- 
ity of which are hosted by disk galaxies, are primarily radio- 
quiet objects. Using high-resolution optical and radio measure- 
ments of a well-defined set of Seyfert 1 galaxies, they showed 
that the nuclear R parameter places more than 60% of the ob- 
jects in the category of radio-loud sources (log/? > 1). The 
present study proceeds in a similar spirit. In order to assess 
the relative radio power, we constructed nuclear measurements 
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of R' , which differs from R only in that the optical continuum 
luminosity was obtained indirectly through the H/3 luminosity. 
When studying nearby galaxies with low-luminosity nuclei, we 
reiterate the importance of using high-resolution data to prop- 
erly define the nuclear R parameter — this is the quantity that 
is relevant for comparison with luminous AGNs and quasars, 
whose nonstellar nuclei dominate the integrated emission. We 
disagree with Laor's (2000) use of single-dish radio measure- 
ments to evaluate the radio loudness of the objects in the sam- 
ple of Ho (1999c). Laor also questioned the utility of the stan- 
dard R parameter as apphed to low-luminosity AGNs because 
he suspected that the optical bolometric correction may be ex- 
ceptionally large for these objects. The characteristic weakness 
of the "big blue bump" in the spectral energy distributions of 
low-luminosity AGNs (Ho 1999c, 2001; Ho et al. 2000) indeed 
does lead to a larger value of Lboi/LB than is typically seen in 
higher luminosity sources, but this difference is only a factor of 
^2 (see § 3.1), which is insufficient to alter the main conclu- 
sions. 

4. We find a striking inverse correlation between R' and 
A = Lboi/LE. Since A varies as a function of the mass accre- 
tion rate M, the most straightforward implication of this result 
is that the relative radio power increases with decreasing M. 
A significant fraction of the strongly active sources have high 
accretion rates (log A ^ -2) and are radio quiet (log/?' < 1), 
whereas nearly all of the weakly active objects are starved for 
fuel (all log A < -2) and are radio loud (log/?' > 1). 

5. The systematic dependence of /? on A and the tendency 
for local galaxies to have nuclei which are both underlumi- 
nous and radio loud are quaUtatively consistent with predic- 
tions from accretion-disk theory. The amount of fuel avail- 
able in the centers of present-day galaxies is plausibly quite 
low. As the accretion rate falls below a critical threshold of 
M « O.Io^Me « O.OIMe, where a is the standard viscosity 
parameter (assumed to have a value ~0.3), the accretion flow 
makes a transition to an optically thin, two-temperature ADAF 
(Narayan et al. 1998b). Under these conditions, the low den- 
sity and low optical depth of the accreting material lead to 
inefficient cooUng, the radiative efficiency is much less than 
the canonical value of 10%, and thus the resulting luminosity 



is low. Moreover, ADAFs produce generically "radio-loud" 
spectral energy distributions, for two reasons. First, cyclo- 
synchrotron emission at radio wavelengths provides energet- 
ically important cooling. And second, the broad-band spec- 
trum, by definition, lacks the big blue bump usually attributed 
to thermal emission from the optically thick, geometrically thin 
disk (Shields 1978; Malkan & Sargent 1982). ADAFs are nat- 
uraUy bright in the radio and dim at optical/ultraviolet wave- 
lengths: both conspire to boost R. The optical/ultraviolet com- 
ponent from an ADAF comes from inverse Compton scatter- 
ing of the cyclo-synchrotron photons, and its strength increases 
sensitively with rising M/M-^n (see e.g., Fig. 1 of Mahade- 
van 1997). This quaUtatively explains the inverse correlation 
between R and A. Furthermore, Rees et al. (1982) have sug- 
gested that the vertically thick structure of the "ion torus" may 
help facilitate the coUimation of radio jets. Finally, we note 
that the majority of the weakly active objects in our sample 
have values of A ^ lO""*- 10"^, comfortably below the thresh- 
old within which ADAFs operate. Although the above general 
arguments need to be confirmed with more quantitative calcula- 
tions, a number of authors have invoked ADAF models to fit the 
spectral energy distributions of some of the objects in our sam- 
ple (Sgr A*, Manmoto, Mineshige, & Kusunose 1997, Narayan 
et al. 1998a; NGC 4258, Lasota et al. 1996, Chary et al. 2000; 
M81, Quataert et al. 1999; M87 and NGC 4649, Di Matteo et 
al. 2000). 
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